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Tema 3: Descripción numérica y gráfica de datos. 





En este tema: 

1.- Estadística básica. 

2.- Distribuciones de Probabilidad. Generación de variables aleatorias. 
3.- Tablas de frecuencias. 

4.- Medidas de localización, dispersión y forma. 

5.- Descripción gráfica de datos en R. 


6.- Gráficos para datos discretos. 
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Tema 3: Descripción numérica y gráfica de datos. 





También: 
7.- Gráficos para datos continuos. 
8.- Representación de datos multivariantes. 
9.- Gráficos para estudiar la distribución de unos datos. Estimación de 
densidades. 
10.- Parámetros gráficos. Cambios permanentes: la función par(). 


11.- Exportando gráficos. Dispositivos gráficos (device drivers). 
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1 Estadistica DÍsica. 


e Los datos obtenidos cuando realizamos cualquier experimento 
presentan variabilidad: 
» el peso de un bebe al nacer varía 


> la cantidad de lluvia recogida en un dia en una determinada zona varía 
> la altura de una planta sometida a dos tipos de abono varía, etc. 


o La Estadística es una disciplina que se ha desarrollado en respuesta a 
los experimentadores cuyos datos exhiben variabilidad. 

e Los conceptos y métodos de la estadística nos permiten describir la 
variabilidad, planificar la investigación teniéndola en cuenta y analizar 
los datos para extraer el máximo de información de los mismos así 
como determinar la fiabilidad de las conclusiones que podamos obtener 
a partir de estos datos. 

o Ya sabemos que R es un lenguaje que permite implementar técnicas 
estadísticas. 

o En este tema vamos a ver cómo podemos utilizar R para analizar los 
bancos de datos que habitualmente se nos presentan al trabajar. 
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L + ESTAdÍstica DÍSICa, 


o Variable —> Característica de interés. 


o Muestra Observada —> Conjunto de valores de la variable observados 
obtenidos de manera homogénea. 
o Tamaño muestral —> Número de datos observados. 
o La manera de describir la muestra (nuestros datos) depende del tipo 
de atributo: 
> Cualitativo —> Intrínsecamente no tiene carácter numérico 
(categórica) 
* Nominal (sin orden entre los valores): Sexo 
* Ordinal (con valores ordenados): Nivel de estudios 
> Cuantitativo —> Intrínsecamente numérico 


* Discreto (cantidad finita o numerable de valores): Número de hijos 
* Continuo (valores en toda la recta real): Altura 
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2.- Distribuciones de Probabilidad. 


o R tiene las distribuciones de probabilidad más comunes implementadas 
en la librería BASE. En otras librerías disponemos de otras tantas. 


o Para cada una de ellas (distrib), disponemos de 4 versiones: 





generador de numeros aleatorios rdistrib 
función densidad /probabilidad ddistrib 
función distribución pdistrib 


función inversa distribución (cuantiles)  gdistrib 





Ejemplo 







x.norm <- rnorm(500) + Simulación de 500 datos normales 
2*pt(-2.43, df = 13) + P-valor de dos colas de una t (13) 
af(0.99, 2, 7) + Percentil 99 de una distribución F(2, 7) 
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Distribuciones de probabilidad en la librería BASE. 





Función 

Normal 
exponencial 
gamma 

Poisson 

Weibull 

Cauchy 

beta 

t de Student 

F (Snedecor) 
Pearson ye 
binomial 
geométrica 
hypergeométrica 
logistica 
lognormal 
binomial negativa 
uniforme 








Utilidad 

rnorm(n, mean=0, sd=1) 

rexp(n, rate=1) 

rgamma(n, shape, scale=1) 
rpois(n, lambda) 

rweibull(n, shape, scale=1) 
rcauchy(n, location=0, scale=1) 
rbeta(n, shapel, shape2) 

rtín, dí) 

ríf(n, df1, df2) 

rchisqí(n, df) 

rbinom(n, size, prob) 

rgeom(n, prob) 

rhyper(nn, m, n, k) 

rlogis(n, location=0, scale=1) 
rlnorm(n, meanlog=0, sdlog=1) 
rnbinom(n, size, prob) 
runif(n, min=0, max=1) 
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3.- Tablas de frecuencias para variables categóricas. 


o Hemos visto un factor es un vector utilizado para especificar una 
clasificación discreta de los elementos de otro vector de igual longitud, 


e y que en R existen dos tipos de factores (variables categóricas): 


> Nominales: No existe jerarquía entre ellos (p.e., colores) 
> Ordinales: Existe jerarquía entre ellos (p.e., grupos de edad) 


o Del mismo modo, dos factores definen una tabla de doble entrada, y 


así sucesivamente. 


o La función table() calcula tablas de frecuencias a partir de factores 


de igual longitud. 


o Si existen k argumentos categóricos, el resultado sería una variable 
k-indexada, que contiene la tabla de frecuencias. 


Ejemplo 
x <- as.factor(1:5) 
table(x) 
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4.- Medidas de localización, dispersión y forma para 
variables cuantitativas continuas. 


o La forma más sencilla de empezar a describir unos datos cuantitativos 
es realizar un resumen estadístico. 


o En la página siguiente tenemos un listado de las más habituales para 
describir la localización y la dispersión. 


o Para analizar la curtosis y la asimetría de unos datos podemos utilizar 
dos funciones de la librería e1071: skewness() y kurtosis(). 


Ejemplo 









x <- rgamma(50,1,3) 

summary(x); fivenum(x) 

mean(x); median(x); quantile(x); quantile(x,c(0.35,0.9)) 
sd(x); var(x); range(x); IOQR(x) 

mas): vic mn). Se ltrca mal pun (es ro l01) 
max(x); which.max(x); x[which.max(x)]; pmax(x[4:8],x[2:6]) 


9/33 
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Medidas de localización y dispersión más habituales. 





Función Utilidad 
sum(..., na.rm=FALSE) Suma 
max(..., na.rm=FALSE) Máximo 
min(..., na.rm=FALSE) Mínimo 


which.min(x) 

which.max(x) 

pmax(...,na.rm=FALSE) 

pmin(...,na.rm=FALSE) 

cumsum(x), cumprod(x) 

cummax(x), cummin(x) 

mean(x, trim=0, na.rm=FALSE) 
weighted.mean(x,w,na.rm=FALSE) 
median(x,na.rm=FALSE) 
quantile(x,prob=(0,0.25,0.5,0.75,1),na.rm=F) 


Posición del máximo 
Posición del minimo 
Máximo en paralelo 
Minimo en paralelo 

Sumas y prods acumulados 
max?s y min?s acumulados 
Media 

Media ponderada 

Mediana 

Cuantiles 





fivenum(x, na.rm=FALSE) 


summary(x, na.rm=FALSE) 


5-Tukey: min, lower-hinge 


mediana, upper-hinge, máximo 
min,1c,mediana,media,3c,max 








IOR(x, na.rm=FALSE) 

range(...,na.rm=FALSE, finite=FALSE) 
var(x, y=x, na.rm=FALSE, use) 

sd(x, na.rm=FALSE) 
mad(x,center,constant=1.4426, na.rm=FALSE) 





Rango inter-cuartilico 
Rango 

Varianza 

Desviación Tipica 


Desviación mediana absoluta 
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Ejemplo 

library(e1071) 

x <- rgamma(50,1,3) 

moment (x,2,center=F) + momento no centrado de orden 2 


Consideramos dos distribuciones asimétricas (Betas) y las vamos a comparar con 
la normal que es simétrica: 

nsim<-5000 

s1<-skewness(rbeta(nsim,2,3)) 

s2<-skewness(rbeta(nsim,3,2)) 

s3<-skewness(rnorm(nsim,0.5,0.5)) 

s1;s2;s3 


Consideramos ahora una distribución normal y una Student, más achatada, y las 
comparamos: 

k1<-kurtosis(rnormí(nsim)) 

k2<-kurtosis(rt(nsim,3)) 

k1;k2 
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5.- Descripción gráfica de datos en R. 


Las gráficas son la mejor forma de simplificar lo complejo. Un buen 
gráfico suele ser más accesible que una tabla. Sin embargo es muy 
importante tener claro qué gráfico queremos hacer. 


Las facilidades gráficas de R constituyen una de las componentes más 
importantes de este lenguaje. 


R incluye muchas y muy variadas funciones para hacer gráficas 
estadísticas estándar: desde gráficos muy simples a figuras de gran 
calidad para incluir en artículos y libros. 


Permite además construir otras nuevas a la medida del usuario 
(aunque a veces hacer cosas simples no es fácil). 


Permite exportar gráficas en distintos formatos: PDF, JPEG, GIF, etc. 
Para ver una demo de gráficos con colores: demo (graphics). 


Aquí únicamente veremos algunas de todas las posibilidades. Más 
detalles en el libro R Graphics de Paul Murrell. 


Otra alternativa: ggplot. 
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Gráficas en R 


o R tiene dos sistemas de producir gráficos: 

» El tradicional, que es el que veremos principalmente 

> Gráficos Trellis (paquete Lattice) del que veremos algunos ejemplos 
o Podemos dividir los comandos para efectuar las gráficas en tres grupos: 


> Funciones para crear gráficas de alto nivel, es decir ya programadas y 
que admiten diferentes posibilidades. 

> Funciones de bajo nivel, que permiten un control más fino del dibujo y 
permiten crear gráficas a medida. 

> Funciones para el uso interactivo, para extraer información de una 
gráfica o una modificación mediante el ratón. 
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Función plot O 


El procedimiento gráfico de alto nivel más habitual para dibujar datos es 
plot (). 


Ejemplo 

x<-(0:65)/10 

y<-sin(x) 

plot (x) 

plot(x, y) 

plot(x, y, main="Función Seno") 
z<-cos(x) 

windows() tHCrea una ventana nueva 
plot(x, z, main="Función Coseno") 
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Opciones de la función plot () 


Algunas de las más útiles 
o main: Cambia el título del gráfico 
sub: Cambia el subtítulo del gráfico 
type: Tipo de gráfico (puntos, líneas, etc.) 
xlab, ylab: Cambia las etiquetas de los ejes 
xlim, ylim: Cambia el rango de valores de los ejes 


lty: Cambia el tipo de línea; 1wd: Cambia el grosor de línea 


col: Color con el que dibuja 





Ejemplo 
plot(x, y, main="Seno", type="1") 

plot(x, z, main="Coseno", lty=2, col="red", type="1") 
plot(x, z, main="Coseno", lty=3, col="blue", type="1", 
xlim=c(0, 2), ylab="cos(x)") 
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Procedimientos de bajo nivel 


Hay una serie de funciones que permiten dibujar sobre una gráfica ya 
creada. 


Los más habituales 







points(x, y, ...): Dibuja una nube de puntos 
lines(x, y, ...): Dibuja una línea que une todos los puntos 
ablines(): Dibuja una línea recta dada la interc. y pendiente 


polygons(x, y, ...): Dibuja un polígono cerrado 


text(x, y, labels, ...): Escribe texto en unas coordenadas 





Ejemplo 
plot(x, y, main="Funciones seno y coseno", type="1") 
lines(x, z, col="blue", lty=2) + col=4 es equivalente 

Lex (Ec (000. 0) 0. 1d). Mlabele=e (ela). cos (a ya 
coléc(blace e "blues )) 
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Leyendas 


Descripción 
La función legend(x, y, legend, ...) permite añadir leyendas a un 
gráfico: 

o x,y : Esquina sup. izda. de la leyenda 

o legend: Texto de la leyenda 

o bty: Tipo de borde (“n” para omitir) 





Ejemplo 
plot(x, y, main="Funciones seno y coseno", type="1") 
lines, 2. colo blue. LLy>2) 

legend(x=3, y=1, legend=c("sin(x)", "cos(x)"), lty=c(1,2), 
colec blecrtbine in) 
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Funciones gráficas interactivas. 


En R existen una serie de funciones que permiten completar los gráficos de 
manera interactiva por parte del usuario 
Descripción 
o identify(x, y, etiquetas) identifica los puntos con el ratón y 
escribe la correspondiente etiqueta. 


o locator() devuelve las coordenadas de los puntos. 





Ejemplo 








plot(x, y, main="Funciones seno y coseno", type="1") 
lines(x, z, col=2, lty=2) 
legend (locator(1),legend=c("sin(x)","cos(x)"),lty=c(1,2),co1=c(1,2)) 


x <- 1:10; y <- sample(1:10) 
nombres <- paste("punto", x, ".", y, sep ="") 
plot(x, y); identify(x, y, labels = nombres) 
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6.- Gráficos para datos discretos. 


Los más habituales 


Para representar variables categóricas o cuantitativas discretas (con pocas 
clases): 


o Diagramas de puntos: dotplot () 


o Diagramas de barras: barplot () 


o Diagramas de quesos: pie() 





Ejemplo 

library (lattice) 

x <- rbinom(100,5,0.3) 

dotplot (table (x),horizontal=F) 
par(mfrow=c(2,2)); plot(x,type="h") + Diagrama de puntos 

barplot (table(x),col=rainbow(length(table(x)))) + Diagrama barras 
pie(table(x)) + Diagrama de quesos 





DET RATE NES UVA) Comp. y Prog. en R: Tema 3 19 / 33 


7.- Gráficos para datos continuos. 


Los más habituales 
o Diagramas de cajas: boxplot () 


o Diagramas de tallo y hojas: stem() 


o Diagramas de puntos: stripchart() 





Ejemplo 















y <- rnorm(100); y.f <- rbinom(100,5,0.3) 

$ Gráfico de tallos y hojas 

stem(y) 

par (mfrow=c(2,2)); m<-mean(par("usr'")[1:2]) + medidas ventana usuario 

+ Diag. de cajas 

boxplot (y); boxplot(y-y.f); boxplot(split(y,y.f),col="cyan") 

$ diagrama de puntos, tres métodos 

stripchart (y); text(m, 1.04, "stripchart método overplot") 

stripchart (y, method="jitter",add=T,at=1.2); text(m,1.35,"método jitter") 
stripchart (round (y, 1) ,method="stack",add=T,at=0.7); text(m,0.85,"método 
stack") 
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8.- Representación de datos multivariantes. 


Cuando queremos representar varias variables conjuntamente para detectar 
relaciones entre ellas, disponemos de diversos tipos de gráficos: 


Los más habituales 
o Gráficos de tendencias para tablas de contingencia: dotchart () 


o Gráficos de dispersión: plot () y pairs() 


o Gráficos condicionados: coplot (). 





Ejemplo 
$ Gráficos de tendencias para tablas de contingencia 
data(VADeaths) 

dotchart(VADeaths, main = "Death Rates in Virginia - 1940") 
H Gráficos condicionados 

data(quakes) 

coplot(lat=1long | depth, data = quakes) 
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Ejemplo 















$ Gráficos de dispersión para revisar relaciones entre variables 
X <- matrix(rnorm(1000), ncol = 2); colnames(X) <- c("a", "b") 
plot (X) 


X <- matrix(rnorm(1000), ncol = 5) 
colnames(X) <- c("a", "id", "edad", "loc", "peso") 
pairs(X) 


data (iris) 

razas<-unclass (irisfSpecies) 

plot(iris[1:2],pch=21,bg = c("red", "green3", "blue") [razas]) 
pairs(iris[1:4], main = "Anderson?s Iris Data - 3 species", pch = 21, bg 
= c("red", "green3", "blue") [razas]) 


data(swiss) 
pairs(swiss, panel = panel.smooth, lwd = 2, cex= 1.5, col="blue") 
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9.- Gráficos para estudiar la distribución de unos datos. 


Cuando queremos estudiar cual es la posible distribución de unos datos 
disponemos de diferentes funciones: 
Los más habituales 

o Histogramas: hist () 

e Gráficos qq: qqaplot(), qqnorm() y aqline(). Dos posibles usos: 


Comparación de cuantiles empíricos versus cuantiles teóricos: para 
comprobar si los datos se parecen a una determinada distribución 
Comparación de dos distribuciones empíricas entre sí 


o Estimación de la función de distribución empírica: ecdf () 


o Estimación kernel de la función de densidad: density() 





Ejemplo 
$ Histogramas 
y<-rnorm(500); hist(y); hist(y,5) 
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Gráficos qq 


Ejemplo 














y<-rnorm(500) 
+ Comparación de los cuantiles muestrales con los de una Normal 
aqnorm(y); qqline(y) 

+ Comparación de los cuantiles muestrales de dos muestras 

y or (0073) 

aqplot(y,y.t,xlab="Dist. Normal", ylab="Dist. St(3)"); qqline(y) 
+ Comparar cuantiles muestrales con los de una distribución dada 
library (lattice) 
agmath(y,distribution=function(p)1qt(p,df=5))) 
agmath(y,distribution=function(p)1qgamma (p,shape=3,rate=5))) 


DET RATE ENEE (UVA) Comp. y Prog. en R: Tema 3 24 / 33 


Función de distribución empírica. 


Ejemplo 

x<-rlogis(500,2,3) 
Fn.x<-ecdf (x); summary(Fn.x) 
plot (Fn.x,main="Función Distribución Empírica") 
plot(Fn.x,add=T,verticals=T,col.v=2,c01.h=4) 
x.seq<-seq(min(x),max(x),length=4) 

$ Para evaluar Fn en cualquier punto 
points(x.seq,Fn.x(x.seq),pch=21,bg="red3") 
abline (v=x.seq,col="red3") 
abline (h=Fn.x(x.seq),col="red3") 
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Estimación de funciones de densidad. 


Ejemplo 
x<-rgamma(500,4,3) 
hist(x,prob=T) + prob=T equivale a freq=F 


$ Estimador kernel de la densidad 
lines(density(x)) 

bw.x<-density(x) $bw 

bw.x $ amplitud de la banda 

$ Podemos modificar la banda para la estimación 
lines (density(x,bw=bw.x/2),co1=2) 

lines (density(x,bw=bw.x*2),co1=4) 
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Representación en 3D. 


Cuando queremos representar una función bivariante disponemos de 
diversos tipos de gráficos: 


Los más habituales 
o Gráficos en tres dimensiones: image () 
e Gráficos de contorno: contour (). Permite añadir líneas de nivel. 


o Las librerías MASS y ks tienen funciones para estimar kernels 
bivariantes. 





Ejemplo 
x<-seq(-1,1,0.05); y<-seq(-1,1,0.05) 
E ss tunctionta, y) cos) /(1 + EN2) 
a ote E) 

image (x,y,z); contour(x,y,z,add=T) 
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10.- Parámetros gráficos. La función par() . 


o Al igual que en la función plot (), podemos controlar casi todos los 
aspectos de una gráfica mediante los parámetros gráficos. 

o R dispone de una lista (unos 70) de parámetros gráficos que controlan, 
por ejemplo, el color, tipo de línea, grosor, tipo de punto, justificación 
del texto, tamaño, etc. 

o Vale la pena leer con calma la ayuda de par. 


o Cada parámetro tiene un nombre (p.e. “col” para color) y un valor. R 
dispone de colores de los más variado. Ver el archivo UsingColorInR. 


o Cada vez que hacemos una gráfica se abre un dispositivo gráfico con 
una lista de parámetros que tiene unos valores iniciales por defecto. 
e Los parámetros gráficos pueden fijarse: 


> de forma permanente para ese dispositivo mediante la función par () 
> O temporalmente en las llamadas a las funciones gráficas incluyéndolos 
en la lista de argumentos (si lo permiten). 
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Funcionamiento de la función par). 


eo Sin argumentos devuelve una lista con los parámetros y sus valores en 
activo. 


Ejemplo 
par () | 


e Con un argumento, vector de caracteres, con los nombres de algunos 
parámetros, devuelve una lista con los parámetros y sus valores en 
activo. 


Ejemplo | 
pario col 


e Con nombres de parámetros = valor, establece los nuevos valores 


Ejemplo 
par(col=4,1ty=2) | 
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Ejemplos de parámetros gráficos 


Colocar varias gráficas en una ventana 
Los siguientes parámetros permiten diseñar el número de gráficas en cada 
dispositivo gráfico 
o mírow: N* de filas y columnas en la ventana. Los huecos se rellenan 
por filas. 


o mícol: Idem pero se rellena por columnas. 





Ejemplo 
x<-(0:65)/10 
y<-sin(x) 
par(mfrow=c(1,2)) + probar también c(2,1) 

plot(x, y, main="Seno", type="1", ylab="sin(x)") 
plot(x, z, main="Coseno", type="1", lty=2, col="red", 
ylab="cos(x)") 
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11.- Dispositivos gráficos (device drivers). 
o El funcionamiento de los diferentes dispositivos gráficos depende 
mucho del entorno de trabajo. 
o Por defecto, cuando realizamos la primera gráfica, R abre un 
dispositivo gráfico. 
o Podemos abrir nuevas ventanas gráficas (p.e. en el entorno Windows 


llamando a la función windows ()). Con ello tendremos varios 
dispositivos donde dibujar. 


o Para cerrar un dispositivo abierto utilizamos dev.off (). Si no 
tenemos claro cual cerrar, la función dev.1list() nos puede ayudar a 
saber qué dispositivos hay abiertos y que numeración tienen. 


o Siempre hay uno activo, podemos saber cuál es con dev.cur(). Si 
queremos activar otro podemos utilizar dev.set (). 


e Con la opción histórico grabando activa R nos permite disponer de 
todos los gráficos e ir accediendo al resto con Av.Pág. y Re.Pág. 


o RStudio permite una gestión muy sencilla de los gráficos realizados. 
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Exportando gráficos. 


o Para guardar una gráfica, podemos copiar y pegar desde la ventana 
gráfica a un tratamiento de textos que los permita. 

o Desde el menú Archivo -> Guardar como podemos guardar la 
gráfica como un fichero metafile, pdf, png, bmp, postscript, tifo jpg. 

e Sin embargo esta opción no es la mejor ya que no tenemos control 
sobre la propia gráfica y como queda guardada. Sobre todo a nivel de 
escala. 

e Lo mejor es enviar directamente la gráfica a un dispositivo (pdf, 
postscript, etc.) utilizando funciones como pdf () o postscript(). 


Ejemplo 
pdf ("prueba.pdf", paper="special", width=13, height=7) 
hist (x<-rnorm(100),prob=T) 

dev.ofíf() 
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Licencia de este material ] 
Más info: http: //creativecommons.org/licenses/by-sa/3.0/es/ 


Usted es libre de: 


copiar, distribuir y comunicar públicamente la obra 
hacer obras derivadas 


Bajo las condiciones siguientes: 












Reconocimiento. Debe reconocer los créditos de la obra de la mane- 


(4) ra especificada por el autor o el licenciador (pero no de una manera 
que sugiera que tiene su apoyo o apoyan el uso que hace de su 









obra). 
Compartir bajo la misma licencia. Si transforma o modifica esta 
obra para crear una obra derivada, sólo puede distribuir la obra 
resultante bajo la misma licencia, una similar o una compatible. 
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